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The backreaction of structure formation influences the cosmological evolution equation for the ho- 
mogenous and isotropic average metric. In a cold dark matter universe this effect leads only to 
small corrections unless a substantial fraction of matter is located in regions where strong gravi- 
tational fields evolve in time. A "cosmic virial theorem" states that the sum of gravitational and 
matter pressure vanishes and therefore relates the average kinetic energy to a suitable average of 
the Newtonian potential. In presence of a scalar "cosmon" field mediating quintessence, however, 
cosmology could be modified if local cosmon fiuctuations grow large. We speculate that this may 
trigger the accelerated expansion of the universe after the formation of structure. 



PACS numbers: 98.80.Cq, 95.35. -fd 
I. INTRODUCTION 

Recent observations of the Hubble diagram for super- 
novae la indicate that the expansion of the universe may 
be accelerating in the present epoch [|]. In this case 
the previous decrease of the Hubble parameter H ^ 
must have slowed down just in the last couple of bil- 
lion years and an|-pbvious question asks: Why does this 
slowdown happerJJ "just now"? A possible explanation 
would be a cosmological constant which sets a mass scale 
A « (10~'^eV)'' and therefore also a corresponding time 
scale ~ (A/M2)-i/2^ ^ith Mp « 10^^ GeV the Planck 
mass. Since it seems to be very hard to understand 
the origin of the tiny mass scale A theoretically, one is 
tempted to look for alternatives. A possible scenario are 
models of quintessence 0, They are based on the 
time evolution of a scalar field ~ the cosmon - which is 
of cosmological relevance today. In the simplest viable 
models, however, the characteristic time for the onset of 
acceleration is put in by hand in-the form of the effective 
scalar potential or kinetic termH. This does not always 
need tremendous fine-tuning of the order of 100 digits as 
for the case of the cosmological constant. Indeed, there 
are models where it is sufficient to tune parameters on the 
level of percent to permille. We feel, nevertheless, that 
these ideas would become much more credible if a nat- 
ural solution of the "why now" problem could be given. 
As one possibility the relevant time scale may be linked 
to a natural small number arising from a fundamental 
theory. A recent proposal in this direction involves the 
properties of a conformal fixed point [Q . As an alterna- 
tive, some event in the more recent cosmological evolu- 
tion could have "set the clock" to trigger the acceleration 
at present. An idea in this direction Q| - "fc-essence" - 
tries to use the transition from a radiation-dominated 
to a matter-dominated universe in order to set the clock. 



^On a logarithmic scale as relevant for cosmology the last 
few billion years are more or less the "present" epoch. 
^For more recent examples see W, [h], H. 
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Here we explore if structure formation could have induced 
a change of the pace of expansion. 

One of the most striking qualitative changes in the re- 
cent history of the universe is the formation of structure. 
For most of the cosmological evolution the universe was 
homogenous to a high degree. Looking at the sky today 
we see, however, strong inhomogeneities in the form of 
stars, galaxies and clusters on length scales sufficiently 
small as compared to the horizon. Could the emergence 
of the inhomogeneities set the clock ||^ for the present ac- 
celeration? In order to answer this question, we have to 
understand how inhomogeneities on the scales of clusters 
or smaller "act back" on the evolution of the homoge- 
nous "average metric" . (For the purpose of this paper 
we consider formally an average over the present hori- 
zon. More accurately, the supernovae results concern an 
average over a volume corresponding to z « 1.) After 
all, the universe is not homogenous at present and the 
Einstein equations determine the metric in presence of 
these inhomogeneities. One can still formulate a type of 
"macroscopic Einstein equation" for the average metric 
which, by definition, can be considered as homogenous. 
The macroscopic equation simply obtains by averaging 
the "microscopic Einstein equation" . In this averaging 
procedure the "backreaction" of the inhomogeneities ap- 
pears in the form of "correction terms" in the macro- 
scopic Einstein equation |lO| . In models of quintessence 
this holds also for the macroscopic evolution equation for 
the scalar field. 

It is the aim of the present paper to estimate the size 
and therefore the relevance of the backreaction effects. 
For this purpose we express in sect. 2 the backreaction 
effects in terms of a "gravitational energy density" pg 
and corresponding "pressure" Pg. It is obvious that pg 
and Pg are relevant only if they are not tiny as compated 
to the energy density p in radiation or matter. A very 
rough estimate shows that in early cosmology the effects 
of Pg and Pg are indeed completely negligible. One the 
other hand, once stars and galaxies have formed the ratio 
Pg/ p is not many orders below one any more, and a more 
detailed investigation becomes necessary. In sect. 3 we 
evaluate pg and Pg in terms of the correlation function 
for the local energy momentum tensor of matter and ra- 
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diation. This form exhibits clearly the relation of these 
quantities to the inhomogeneities. 

In sect. 4 we attempt a quantitative estimate 
for a standard cold dark matter universe (without 
quintessence, but possibly in presence of a cosmological 
constant). We find that the effects of inhomogeneities on 
the scales of stars and galaxies are small, they contribute 
typically Pg/p ~ 10~^. A typical contribution from in- 
homogeneities on the scales of clusters is pg/p ~ 10"*. 
These estimates hold, however, only if the fraction of 
matter in regions of strong gravitational fields like black 
holes or the center of galaxies is small. In sect. 5 we ad- 
dress the backreaction effects from black holes and similar 
objects. The gravitational energy density pg can indeed 
be large. Nevertheless, the combined energy momentum 
tensor for gravitational and matter contributions behaves 
as for a nonrelativistic gas if the objects are static. We 
conclude that for a cold dark matter universe the back- 
reaction effect could play a significative role only if a 
substantial fraction of matter is found in regions where 
strong gravitational fields evolve in time. This does not 
seem to be very likely. 

In models with quintessence the situation could change 
dramatically, but only if the inhomogeneities in the cos- 
mon field are substantial. The "gravitational backreac- 
tion" Pg,Pg is then supplemented by a "cosmon backre- 
action" PcPc due to the cosmon fluctuations. We dis- 
cuss a simple collection of static and isotropic cosmon 
lumps in sect. 6. This would behave similar to black 
holes. We argue that for more general, in particular non- 
static, cosmon fluctuations the fine cancellation between 
Pc and Pg observed in the static isotropic solutions may 
not be maintained. In particular, it seems conceivable 
for a "cosmon dark matter" scenario [pl| that Pc/ p and 
for Pg/p are of order unity. Under this condition it would 
become quite likely that the formation of structure would 
lead to a qualitative change in the evolution equation for 
the average metric. One would expect deviations from 
H ^ t^^ once structure has formed. We summarize our 
conclusions in sect. 7. 



II. THE INFLUENCE OF STRUCTURE 
ON THE COSMOLOGICAL EQUATIONS 

After the formation of structure the universe does not 
remain homogenous on small scales. Nevertheless, we 
believe that homogeneity and isotropy are realized on 
large scales and describe the cosmological evolution by 
a Robertson- Walker metric. The true metric g^jy of the 
universe has to reflect the inhomogeneities due to stars, 
galaxies and clusters. Therefore the homogenous cosmo- 
logical metric can at best be interpreted as some type of 



average metrical g^i/ =< g^n >. This situation introduces 
a mismatch in the standard treatment of the cosmologi- 
cal Einstein equations. On the right-hand side one uses 
the average of the energy momentum < >, whereas 
for the left-hand side one employs the Einstein tensor 
formed from the average metric, g^i^. The correct aver- 
aged Einstein equation involveso, however, the averaged 
value of the Einstein tensor 

1 1 

< Rf.. - ^Rg^c. >- ^ < t^, > (1) 

The difference between the averaged Einstein tensor and 
the Einstein tensor formed from the average metric, i.e. 
Rfj-f ~ ^Rg^vi introduces a correction term in the cosmo- 
logical equation for the average metric 

- \r9^.u = ^T^. = V > +TP (2) 

Here the "gravitational correction" to the energy momen- 
tum tensor 

= -2M2 < SG^, > (3) 



SGfi„ = R^u — -^Rgt^u — {Rfii^ — ^Rgfiv) (4) 

reflects the influence of the inhomogeneities. It accounts 
for the "backreaction" of structure formation on the evo- 
lution of the homogenous "background metric" (jfj,^. Ho- 
mogeneity and isotropy of all averaged quantities imply 
that the only nonvanishing components of T^^^ are given 

by 

Tio ^Pa = -2Af2 < 5Goo > 
^^Pg-g,,^-2M^ <6G,,> (5) 

The cosmological equation therefore preserves its form, 
but Too is not given solely by the average of the energy 
density in matter and radiation. It also contains a gravi- 
tational contribution which reflects the imprint of struc- 
ture formation on inhomogeneities of the metric. We 
observe that for a flat background metric g^i, = rj^^ the 
quantity T^^^ represents precisely the definition of the 
gravitational energy momentum densities [ p^ . Our set- 
ting is therefore a straightforward generalization to cos- 
mology. 

At this point some comments about our averaging pro- 
cedure seem in order. Assume that in a given suitable 
gauge - we will later specify a particular one - the detailed 



^The averaging is done here with respect to the background 
metric. See ref. ^] for a recent discussion of averaging pro- 
cedures in a more general context. 

■'We use signature (-,-|-,-|-,-|-), R^^p^ = —dpT^p + ... and 
= MpV(167r) = l/(167rG]v). 
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inhomogeneous geometry of the universe is described by 
the microscopic metric g^i, (x, t) . It is related to the mi- 
croscopic energy momentum tensor t^i, (x, t) by the mi- 
croscopic Einstein equation. 



2R'P 



(6) 



We next take a reference metric 5^,^ of the homoge- 
neous and isotropic Robertson- Walker form ds^ = —dt^ + 
a?(t)dSF. The microscopic metric can be written as 
giiu = g^u + hp^u- Our reference metric defines surfaces of 
fixed t. At any time t we define the average by 



{A) {a 



.t) = ^ J^d^yAiy 



X, t) 



(7) 



with V a very large volume, (typically of horizon size) in 
the comoving coordinates H. We can then fixJ:he scale 
factor a[t) in g^^ self-consistently by requiring tl {h^j^u) ~ 
0. This finally identifies g^i, and (g^u). 

Spatial averaging at fixed t has been proposed by Fu- 
tamase jlj] . We are free to use such an averaged descrip- 
tion - the only physical assumption in our paper concerns 
averaged homogeneity and isotropy of the real universe, 
namely in our paper (h^,^) = and (i^j/) depending only 
on t. A much more subtle point is the question to what 
extent a real observer actually observes the "spatially av- 
eraged" quantities in the way introduced here. Detailed 
studies conclude [|l5| that this may actually be the case 
- we will not address this topic in the present paper. 
For a particular picture of the cold dark matter scenario 
Futamase concludes that backreaction effects are small 
whereas Buchert et al. speculate ||l2| that the influence 
on the cosmological evolution could be substantial, nev- 
ertheless. 

Let us next discuss the general structure of pg and Pg 
(eq. As long as gravity remains weak, one can expand 
in the small inhomogeneities of the metric 



such that 



(8) 
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Here the differential operators D and E involve two 
derivatives acting on /i or They will be computed more 
explicitly in sect. 3. From < hap >= one concludes 
that Pg and pg are quadratic in h, 



Pa = 

Pg = 



~2M^ < E^I^^'Kph^s > 
< Elf'^^haph-fS > 



(10) 



Thus Pg and Pg involve the correlation function for the 
metric and do, in general, not vanish. 

The local variation of the metric reflects the local vari- 
ations of the energy momentum ten^r according to the 
"microscopic" Einstein equation (^) Within the linear 
approximation to eq. (H), namely 
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the metric fluctuations h^j^y are linear in the fluctuations 
of the energy momentum tensor 5tp,j. In consequence, pg 
can also be viewed as the effect of a nonvanishing corre- 
lation function for the fluctuations of the energy density. 
This correlation function can be observed as a galaxy - 
or cluster - correlation function on appropriate length 
scales. In particular, we know that on small scales the 
universe is far from homogenous. As an example, on 
the length scale of the size of stars very dense regions 
(stars) contrast with an almost empty environment. This 
is equivalent to a huge correlation function and brings us 
back to the question: Can the formation of stars or galax- 
ies influence the evolution of the universe as a whole? j-. 
In order to get a first rough estimate of the magnitude^ 
of this "backreaction of structure formation" let us as- 
sume for a moment that the universe consists of ran- 
domly distributcdEI "stars" with radius L or volume 
vl — and density p^. Consider our horizon vol- 

ume V with Ny stars. Since Tqq = p is of the same 
order as < too >— NvvlPl/V and, on the other side, 
< ^00 >= ^vvlp\/V, one has 



<5p^ > _< (too - P? > 



1_ r _ NyVL _ _p_ 

/ ' ^ " V ^ PL 



(12) 



The fraction of volume occupied by stars, /, is indeed a 
tiny number and one concludes that the relative density 
fluctuations are huge! On the other hand, the weak grav- 
itational coupling enters such that the relative size of pg 
as compared to p could still be small. Since the operator 
D in eq. ( pi] ) contains two derivatives a rough estimate 
assumes 



^ Since we average at fixed t it does actually not matter if we 
average over coordinates a{t)x or x if the physical averaging 
volume grows ^ a? . 

®In order to achieve this task we may use the freedom of se- 
lecting a suitable gauge. Of course, {h^u) = is only possible 
if the universe is really homogeneous and isotropic in average. 
In particular, (/i/j^) should not depend on x (after using the 
gauge freedom). 



''"Microscopic" means here the scales of stars or galaxies...! 
These scales are to be compared with the "macroscopic" scale 
of the order of the horizon. 

*Note that this estimate does not account for the total back- 
reaction. Since we want to study here the effects of struc- 
ture formation we can concentrate on a typical wavelength 
well within the horizon. Discussions of the backreaction from 
modes outside the horizon can be found in M, pj. 

^ Stars may be replaced by galaxies or other extended 
objects. 
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M4 



<Sp^ > 



(13) 



and, with a similar dimension argument, 

The critical quantity for the relevance of the backreaction 
is therefore given by the ratio 



(14) 



R 



Pg_ 
P 



AP 



PL 



2/3 1/3 

H Pl 
A/2 



LAP 
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with rriL the mass of the stars. It is suppressed by two 
powers of the Planck mass as expected for a gravitational 
fluctuation effect. On the other hand, the mass m^, of 
the star and its size L are huge in microphysical units. 
Inserting values typical for the sun, = 2 • 10'^^ <? = 
1.1 • 10^^ GeV, L = 7 • lOSi = 3.5 • 10^^ GeV-^ and using 
M = 1.72 • 10^* GeV one finds for main sequence stars 
rriL 



LM2 



10" 



(16) 



Another estimate relates the gravitational backreac- 
tion effect to typical values of the Newtonian gravita- 
tional potential (j) — —hoo/2 in extended objects. In- 
deed, we note that R is proportional to the gravita- 
tional potential at the surface of the star, m^G/L, with 
= 167rM2. Its value for the sun IS 

- = — ^ = 2.12 • 10"*^ (17) 
L 

Similarly, for idealized neutron stars with mass at the 
Oppenheimer-Volkoff limit, rriL = lA-lO'^^g, L = 9.6km 
one has 

These first estimates are, perhaps surprisingly, not much 
below one (as could have been expected from the fac- 
tor M-2). A more detailed investigation, including the 
various proportionality constants and the distribution of 
objects with different L and M^, becomes necessary. 

Before doing so, it is instructive to discuss a few quali- 
tative aspects of the dependence of the ratio R on L and 
PL- (1) The ratio Pg/p is independent of p. It therefore 
shows essentially no time-dependence once the objects 
have condensed with a stationary density and size. (2) 
For a fixed density pL the contribution from smaller ob- 
jects vanishes rapidly. For example, the condensation to 
dust particles or planets is many orders of magnitude too 
small to be relevant. (3) Microphysical objects like nuclei 
play no role for pg (i.e. R « lO^'^^ for a gas of nuclei). In 
early cosmology the contribution of pg is therefore com- 
pletely negligible. (4) For an (elliptical) galaxy consisting 
of vg roughly uniformly distributed stars within a radius 
Lg the density scales pG ~ vg{L/ Lq)^ Pl- (There may 
be some moderate enhancement from dark matter.) For 
a uniform mass distribution in a galaxy the combination 
LqPG = '^g{L/Lg)L'^Pl is changed by a factor vgL/Lg 
as compared to stars. 



III. GRAVITATIONAL ENERGY MOMENTUM 
TENSOR IN COSMOLOGY 

We next turn to a more quantitative discussion of eqs. 
(^ and (p^ ) for the metric inhomogeneities. Since the rel- 
evant length scale for the dominant fluctuations is much 
smaller than the horizon, we can neglect derivatives act- 
ing on (jp^i, as compared to those acting on This 
yields the microscopic field equation up to quadratic or- 
der in the metric fluctuations 

j^\haPSQ Q h J^d Q h -Q Q h 
2 X^fi^i^'^ap ^ ^a^p'^^u ^p^p'^au 



{dph^pd^h^p + 2d^h^pdphp, - dPh^dphp, 

+2d"hP^dchp, - 2d''hp^dPh^, ~ 2dch''Pd^Kp 

~2d^h^pd,h^p + dPh^^d^Kp + dph\d,h^p} 

[id'^hppdM' + Ado,h'^Pdph% 
-d°'hPpdc,h^p - 4.dah"'^dphPp - 2dahP'^dph°'p^ 



2AI 



2 ^^f-LU 



(19) 



Here the indices are raised and lowered with the homoge- 
nous background metric g^j/. The average < SGfj,i, > 
concerns only the part quadratic in /i^j/, since < h^^^ > 
vanishes by definition. It is homogenous (it involves a 
volume integral) and we can therefore perform integra- 
tion by parts for the space derivatives. On time scales of 
the order of the characteristic length scales of the fluctu- 
ations < SGfj^i, > is also essentially static. This allows us 
to perform integration by parts for the time derivatives 
as well, and we obtain 

< SG^, >=l< h-^'^pdji^p + -ih-^^dPOph^, 

~2h"Pdpda,hp, - 2hP^d"da,hp, + 2hppdPda,h% 

-h'^^d^dphP, - h'^^d^dphP^ > 

-^g^. < h^Pd^dpKp + h'^^dPdph^ 

~2h''Pdpdf}hP^ > (20) 

On the other hand, the linear part of the field equation 
relates the metric perturbations to the perturbations in 

d^hpij + dpduh- d^dphP^ - d^dphP^ 
= -^i^*t^--yi'p9p.) 
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(21) 



where we use — d^dp and h = h^^. 

Eqs. (^0|) and (|l|) simplify considerably in the har- 
monic gauge which we adopt from now on 



The linear field equation becomes 



(22) 



1 



1 



and the quadratic metric fluctuations read 



+2hd^h^, > -Ig,,, < h'^Pd^h^p + Ihd^h > (24) 
o 2 

Neglecting graviational waves eq. (|2^ ) has the retarded 
solution 

M^-,r)^^/d3.- ^--^-(-:^-J5--^'') (25) 



47rM2 



|x — x'l 



where r obeys (ir = dt/a. We recover Newton's law for 
the graviational potential (j) = ~^oo/2 for static point 
sources. 

For a distribution of star-like objects the time deriva- 
tives of h^i, involve the peculiar comoving velocities of 
these objects. Since the peculiar velocities are small 
as compared to the speed of light, we can neglect the 
time derivatives in eq. (^ as compared to the space 
derivatives. Furthermore, by virtue of rotation symme- 
try the expectation values involving only one derivative 
in a given space direction vanish and one infers 

Pg ^ -2A/2 < ,5Goo -2Af2 < ^h^PAhap 

8 

+ ^hAh - ^hP^Ahpo + i/iA/ioo > 

2 hp 1 
Pa = ^ < SGI >- 2Af2 < -h'^PAh^p 

+ ^hAh + hp'Ahp, - IhAhl > (26) 
48 6 6 

with A = g^^didj. Sums over all double indices are im- 
plied, with latin indices running from 1 to 3. To leading 
order we only need to keep Stoo such that Ssqq = ^Stoo = 

\5p,5sij = \5pgij and ft-y = /looffy , ^oi = 0, /i = 2/ioo. 
This results in 



Pg = --Af2 < /looA/ioo > 



Pg = < /looA/loo > 



(27) 



and we infer the "equation of state" for the gravitational 
energy momentum tensors of starlike objects 



(28) 



Using ( p5|) we can also express pg in terms of the 

correlation function for the energy density fluctuations 
(y = a3 J ^3^) 

^ d^x [ d^y-J-—5p{x)5p{y) (29) 



327rM2y 



\x-y\ 



(The "retardation" in eq. (|25| ) can be neglected since it 
involves again the peculiar velocities.) It is instructive to 
employ a comoving Fourier basis 



where 



< Sp{k)Sp{k') >= G{k){2n)^S{k + k') 



(30) 



(31) 



The two-point density correlation function G{k) depends 
only on the invariant k'^ = k^. One flnds 



9a' 
SAP 



d^k 
(2^ 



k-^G{k) 



(32) 



From p{~k) = p*{k) one infers that G(fc) is a real posi- 
tive quantity. This implies that pg is positive whereas Pg 
is negative. For small k or long distances G{k) decreases 
rapidly and the fc-integral is infrared finite. The inter- 
esting part comes from large k, where condensed objects 
like stars contribute. On these scales it is convenient to 
switch to physical momenta q = k/a such that 



Pg 
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8Af2 



(2^ 



q-'G(q) 



(33) 



Here we employ 5p{x) — J -^^e^"-'^'-' 5p{q) and < 

5p{q)6p{(f) >— G{q){2Tr)^6{q—q') is the correlation func- 
tion as a function of physical (not comoving) momenta. 
For a given static G{q) the gravitational incoherent en- 
ergy density pg would not depend on the scale factor. 
However, the condensed objects are diluted by the cos- 



mological expansion, and G(g) 



-3 



implies pg 
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similar to the energy density in dark or baryonic matter. 
We conclude that pg is essentially a fixed fraction of the 
energy density of matter : R — pg/ p > Q is independent 
of time. 

Together with the gravitational equation of state (|2^ ) 
this can actually be used for an estimate of corrections 
to the equation of state of matter, p = WmP- Indeed, 
the matter and gravitational energy momentum tensor 
are not separately conserved. Gravitational potentials 
lead to pecular velocities and therefore to nonzero p. In 
other words, the equation of state Wm = holds only 
for "free particles" (ideal dust), i. e. if the gravitational 
interactions are neglected. As for all interacting systems 
we expect corrections. Conservation of the total energy 
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momentum tensor is, of course, exact and implies for 
dR/dt = 



p{l + R) + ZH{p + p + pg+pg) = Q 



(34) 



If p is dominated by massive objects or massive nonrel- 
ativistic particles we can approximate p = pM + 3p/2 
where pM is the contribution of the particle masses. As- 
suming that no masses are added or changed during the 
relevant period in the cosmological evolution we infer 
pM ~ 0'~^iPm = ^'iHpM- Furthermore, if p/p is ap- 
proximately constant this extends to p = SHp. Eq. 
(B4) therefore yields the simple relation 



p + Pg =^0 



(35) 



and we infer an estimate for the pressure of matter which 
is due to the gravitational interactions 



p = WrnP ■ 



Pa u- ^ - 
Pa 27 



(36) 



Since R is small (see the next section for an estimate) 
this amounts only to a tiny correction, justifying our ne- 
glection of pecular velocities. We note that the estimate 
( p5| ) ( psf ) plays the role of a "cosmic viriaLtheorem" since 
it relates the average kinetic energy (p)ll2l to the average 
gravitational potential {pg). 

Turning our argument around we emphasize that the 
relation ( |35| ) implies the "cold dark matter expansion 
law" p ~ a^^ provided R = Q. If we would neglect Pg, a 
nonzero pressure p would correct the expansion accord- 



ing to p ■ 



-3(l-l-!ii„) 



This correction is cancelled by the 



presence of the "gravitational pressure" Pg. Thus backre- 
action effects play a role for the evolution - its role being 
to ensure the validity of the cold dark matter expansion 
law even in presence of pecular velocities or nonzero p. 
Corrections arise only for periods where R or Wm do not 
vanish. 



For small /iqo this "weighs" the energy contrast 5p with 
the Newtonian potential 
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P. = -2 < 



> 



(38) 



For compact objects 5p is almost equal to the local value 
of p. For starlike extended objects the size of their 
own gravitational potential is maximal at the surface, 
4>max — —mG/L. For small (pmax the contribution of 
isolated stars to pg is therefore suppressed by a factor 
~ 10"® as compared to their contribution to p, in accor- 
dance with eqs. (^, ([l^), (|T^). We need, however, also 
the contribution of other stars to <f>. This becomes par- 
ticularly simple in the language ( |37| ) or (38). As long as 
gravity remains weak, we only have to fold any mass con- 
centration with the gravitational potential at the same 
location. Incidentally, this shows that our previous asso- 
ciation of the relevant ratio R = pg/ p with the Newtonian 
potential can be made quantitative 



R 
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< >= 



H<t>5p) 
2p 



(39) 



where ^ S> means an appropriately weighted value of 
(j). This also yields a quantitative value for the pressure 
of matter (and therefore the kinetic energy or pecular 
velocities) according to the "cosmic viral theorem" 



p 



< > 



(40) 



Note that equilibration has not been invoked for this es- 
timate. The cosmic virial theorem follows directly from 
Pg/p =<C ^ /6, 7? = and p/p ~ const.. The average 
kinetic energy density {T)/V = Zp/2 = - < > /4 
may be compared with a virialized gravitationally bound 
system where {T)/V ^ -(0>/2. 

For cold dark matter galaxies the value of the galactic 
gravitational potential in the outer regions, in particular 
the halo, can be estimated from the rotation velocities 



IV. DO STARS AND GALAXIES MODIFY THE 
EXPANSION OF THE UNIVERSE? 

In this section we estimate the size of the backreaction 
effect quantitatively for a standard cold dark matter uni- 
verse. We can use eq. (23) in order to express A/ioo in 
terms of 5p and obtain from eq. (p7|) 



P9 = ^< f^OoOp > 



(37) 



^"The value of pg depends on the precise definition of this 
quantity, e. g. Pg = 2M^ {Go.g"") = Pg - 2M^ {C'i^h'"') = 
(5/9)pg. This does not affect pg and the relation p + Pg = 
whereas the ratios Pg/pg and Pg/p get modified. 



(41) 



Within the halo (r < r/f) the dependence of on r is 
approximately logarithmic 



.In 



(42) 



With Vrot = 0(10""^) we conclude that the galactic poten- 
tial is of similar size as the local potential on the surface 
of a typical star (|l7|). 

Clusters of galaxies, however, have a deeper potential 
well. A typical value for a cluster is 



= -10" 



(43) 



If most matter is found within clusters, this gives an ap- 
proximate lower bound for the gravitational energy den- 
sity 
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Pg^llUP (44) 

We observe that this effect resuhs from the mutual co- 
herent correlations between all the stars in a cluster. The 
dominant length scale of this contribution to the corre- 
lation function (^) is related to the size of the cluster. 

There may still be sizeable contributions arising from 
correlations on smaller scales. The center of the galaxy 
typically contains a region with large gravitational field. 
In this region, however, our linearized analysis does not 
apply any more. A similar statement holds for individual 
black holes outside the center of the galaxy. The precise 
evaluation of these contributions to pg needs a nonlin- 
ear analysis and depends crucially on the question how 
much of the matter in the universe is found in regions 
with a strong gravitational field. We will briefly turn to 
this question in the next section. Only if such "strong 
field contributions" are substantially above the bound 
(^), the gravitational energy density could be relevant 
for the evolution of the universe. On the other hand, 
for a moderate contribution from strong field regions the 
backreaction effect remains small for conventional dark 
matter galaxies and clusters. A value pg/p ~ 10^'^ seems 
to be too small to substantially modify the evolution of 
the universe after structure formation. 



V. CONTRIBUTION OF BLACK HOLES 

For black holes and other regions with strong gravi- 
tational fields the linear analysis of the preceeding sec- 
tions does not remain valid. For an individual black hole 
- or any other static and isotropic object - in a flat 
space-time background the sum of matter and gravita- 
tional energy density is fixed, however, by a conserva- 
tion law. This also holds for the pressure. These laws 
can be expressed in terms of linearized gravity and 
are the analogue of charge conservation in electromag- 
netism. We parametrize the static and isotropic metric 
outside a mass concentration in "isotropic coordinates" 
as ds^ = —B(u)dt^ + C{u)dxdx with = xx. The sum 
of the energy densities and the total pressure are related 
to the functions B and C by the linearized Einstein equa- 
tions, with C — dC /du, etc., as 

p{u) + pg{u) = -2M\C" + -C) 

u 

9 /lf2 2 2 

p{u) + pg{u) = —{C" + -C + B" + -B') (45) 

O IX u 

Also, using Gauss' law, one finds for the integrals over a 
volume with u' < u: 

m{u) =4:TT J du'u'^{p{u')+ pg{u')) = -~8ttA4^u^C' (u) 
P{u) = 47ry" du'u''^{p{u')+pg{u')) 

= \m{u) - ^M^u^B'iu) (46) 



For the Schwarzschild metric the functions B{u) and 
C(u) are given by {G'^ = 16nM^) 

Ciu).{,,^y (47) 

This yields, in particular, m(u — > oo) — m with m the 
total mass of the object related to the Schwarzschild ra- 
dius i?5 = m/(87rM^). Similarly, we observe that the 
integrated pressure vanishes, P{u oo) = 0. 

On a length scale which is large as compared to the 
characteristic size of the objects a collection of static 
isotropic objects - including black holes - can be viewed 
as a collection of point particles with masses mi . The to- 
tal energy momentum tensor T^^i, in eq. (j^ averages both 
the matter and gravitational contributions. If the objects 
are sufficiently distant from each other, this amounts to 
summing mi{u — *■ oo) and Pi (u oo). A collection of 
static isotropic objects behaveSp-therefore like a nonrel- 
ativistic gas with zero pressuralil. In particular, black 
holes that have already formed before structure forma- 
tion - this is the meaning of "static" in a cosmologi- 
cal context ~ behave just as a contribution to cold dark 
matter. Irrespective of the fact that their gravitational 
energy density pg can be substantial, the backreaction ef- 
fect from condensed black holes during or after structure 
formation would not lead to a deviation from the usual 
equation of state. 

The only loophole in the argument that backreaction 
effects can be neglected in a cold dark matter universe 
remains the hypothesis of a substantial contribution from 
black holes forming during or after structure formation. 
For such objects we cannot use the static approximation 
(^). At this stage we cannot exclude that a nonzero, 
perhaps even negative, pressure could play a role for non- 
static regions with strong gravitational fields. 

For a cold dark matter universe we conclude that a 
sizeable influence of backreaction effects is only possible 
if a substantial fraction of the energy density is found 
in regions of strong gravitational fields which evolve in 
time and cannot be described effectively as static ob- 
jects. For being important at the time relevant to the 
supernovae Hubble diagram such a hypothetical evolu- 
tion would have to persist at a redshift z w 1. Discard- 
ing this - perhaps rather unlikely - possibility we find no 
relevant backreaction effect from structure formation in 
a cold dark matter universe. 



^^The vanishing of the pressure including the gravitational 
contribution pg is actually even better obeyed as if pg had 
been neglected. Note that < ffi^ > -|-T®^ is covariantly con- 
served with respect to the background metric p^i^ by virtue 
of eq. (Q) . 
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VI. COSMON FLUCTUATIONS 

Recent cosmological observations suggest the presence 
of a homogenous dark energy component. It has been 
proposed that the dark energy density is time-dependent 
and can be described by the dynamics of a scalar field, 
the cosmon Q . If this "quintessence" scenario ||^ , |^ is 
true, one may also suspect that inhomogeneities in the 
cosmon field could be associated with extended struc- 
tures |l^] . In this section we argue that the backreaction 
effect of structure formation is much stronger in a "cos- 
mon dark matter universe" 1 1 f | than in the standard cold 
dark matter universe. One main reason is the direct con- 
tribution of cosmon fluctuations to the averaged energy 
momentum tensor. One also observes large space-like 
components of the gravitational field in cosmon lumps, 
contributing to large pg and Pg. We underline that the 
material of this section is only relevant if the present local 
fluctuations of the cosmon field are really substantial - a 
possibility that remains speculative as long as no consis- 
tent picture of a cosmon dark matter universe has been 
developed. If the scalar field mediating quintessence re- 
mains homogenous to a high degree in the present epoch, 
its "backreaction" effects are small and can be neglected. 

There are three new ingredients for the backreaction 
in presence of an inhomogenous scalar field: 

(1) Local fluctuations of the scalar field around its ho- 
mogenous background value induce a new contribution 
to the total energy momentum tensor (^. The scalar 
contribution to the local energy momentum tensor 

C = -V{^)g^,, + d^^d,^ - ^^dP^dp^g^, (48) 

yields, after averaging in eq. (|l|) or (^), both a contribu- 
tion from the homogenous background field Cp (the aver- 
age of Lp) and from the inhomogenous local fluctuations 
of the cosmon field Sip = if — (p, namely 



T, 



Here stands for the time variable dark energ; 
mogenous quintessence and corresponds to eq. (M 
Lp replaced by (p 



(49) 

or ho- 
BT) with 



T^o = Ph = Vm + \cp^ 



Tij = Phgz] , Ph 



1 .2 



(50) 



The difference T,'; 
fluctuations (similar to eq. 
in the form 



> 



-T^^ is due to the cosmon 
and can again be written 



00 



Pc 



- Pcgij 



(51) 



It has been discussed in . 

(2) The evolution equation for the background scalar 
fleld also obtains a contribution q''^ from backreaction 
effects llOll [HI 



^ dV _ 
ip + 3Hip+ —{ip) 



(52) 



For cosmon dark matter the "incoherence force" g"^ has 
been discussed in [Q. We note that q'^ can also receive 
a contribution if the cosmon couples to "standard" cold 
dark matter [ p^ . Such a contribution would not be af- 
fected by structure formation. 

(3) The gravitational energy density pg and pressure 
Pg can be enhanced as compared to cold dark matter. 
We discuss this possible effect in a simple model of a 
collection of cosmon lumps ||l^ . We expect that the most 
important features could be present also for more general 
nonlinear cosmon field configurations beyond the specific 
model considered here. 

Let us consider a collection of cosmon lumps (some-of 
them could be associated to some of the galaxiesEJ) 
which can be described in comoving coordinates as 



ip : 



5oo 



9i3 



-{l + J2{BiM--l)} 
a^m,{l + Y,iCdui)-l)} 



9oi = 



(53) 



Here xi is the comoving coordinate of the center of the 
lump £ 



(54) 



and (p{t) is the cosmological background value of the cos- 
mon field ip which leads to homogenous quintessence. 
We assume that the lumps are well separated such that 
6pi,Bi and C'i can be determined from the coupled 
gravity-scalar field equation for a single (spherically sym- 
metric) lump. This system has been discussed in and 
we concentrate on the "halo region" which may give an 
important contribution to the enerrai, momentum tensor. 



In this region we can approximateEj {Rj^ g/e < uj < 



Ci 



B, 



1 



In" 



eui 



R 



Ha 



1 



■In 



Sipi — jeM In Bi 



■In 



In" 



Rhj 



(55) 



^^For our own galaxy a large "cosmon halo" seems unlikely 
in view of the strong distortion of light trajectories jl^ . 

^■^Note that the singularities at ue = RH,t/e correspond to 
pointlike singularities at = in Schwarzschild coordinates. 
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where the scale Rh i can be associated with the radius 



of the halo and 1/|7^| 



is associated to the ro- 



tation velocity of objects in circular orbits within the 
halo. The spherically symmetric solution of the coupled 
gravity-scalar system in empty space has indeed two in- 
tegration constants {Rhj, 7f)- The total mass of the ob- 
ject can be expressed in terms of these constants [0. We 
consider here small values of l/|7f| which correspond to 
realistic rotation velocities and halo extensions of galax- 
ies [0. While Be is close to one except for the vicinity 
of the singularity , we observe that Cg deviates substan- 
tially from one within the halo. As a consequence, hij /a^ 
is of the order one within the halo region and we expect 
substantial contributions to the backreaction effects! 

For cosmon lumps the spacelike components of the en- 
ergy momentum tensor are important. This contrasts 
with stars. For a static lump the time derivative of the 
scalar field vanishes, and one finds for a single cosmon 
lump 



(56) 



These relations are easily generalized to a collection of 
well separated lumps. The cosmon part of the energy 
momentum tensor obeys the equation of state 



Pc 



WcPc 



|(l + 2< AF>/pe) (57) 



Here we have kept only the contribution from the inho- 
mogenous fluctuations and AV is the difference between 
the local value of the cosmon potential and the homoge- 
nous cosmological valualj. For simplicity we will later 
concentrate on the case where < AV > can be neglected 
such that Wc — —1/3. Within the halo the potential con- 
tribution is indeed small. More generally, the background 
potential V{if) is small as compared to the local energy 
densities such that effectively AV > 0. This implies that 
static lumps lead to a negative cosmon equation of state, 
Wc < -1/3. 

We next turn to the gravitational contribution. For a 
spherically symmetric static lump we can again use the 
relations ( ^ ) and (|4^) for the total energy density and 
pressure. In particular, far away from the lump the solu- 
tion approaches the standard Schwarzschild solution p^ ] 
and we infer that the total integrated pressure vanishes 
(uc = Rh/s) 



P{u -> oo) = / dV{pc +Pg+ Pm) = 



(58) 



This implies a cancellation between a negative cosmon 
and positive gravitational contribution! (We also have 
included possible matter {pm,pm) in the "bulk" of the 
object.) A partial cancellation also happens for the en- 
ergy density 



m(u 



) = / dV{pc + Pg + Pm) = m 



(59) 



Indeed, if 7 is large, the Schwarzschild radius Rg — 
TO/(87rM^) « 2Rh/\'^\ is small compared to the halo 
radius Rh and this is equivalent to a substantial cancel- 
lation between pc and pg. 

For an understanding of these cancellations in more 
detail it is instructive to study the gravitational contri- 
bution in the linear approximation. Since the coordinates 
used for the metric ( |53| ) are notiiarmonic, the formulae 
(^) and ( p6|) receive corrections.EJ Being only interested 
in the qualitative features we neglect these corrections 
here. This yields for the gravitational energy density and 
pressure 

Pg = < Ih'^'Ssf,^ + Ihg^'^Ssf,^ + h^'^Ss^o + hSsoo > 
4 o 

Pg = -< ^h^'"Ss^, + ^^hg^'-ds^. 



1 



1 



(60) 



If the potential term can be neglected, the cosmon lumps 
obey (5s^o — such that 

Pg = \< {'2h'^ + hf^) 5s,, >, 

Pa^\<i-^4^'' +hr^5^^,> (61) 



We evaluate the above expression in coordinates 
adapted to the present cosmological time with a — 
1 , g^^ = 5^^ and assume first that only the halo region 
of the lumps contributes effectively to Pg,Pg,Pc a-nd Pc- 
We can therefore evaluate the ratios Pg/pg an pnjpc for 
a single cosmon lump. In the coordinate system (|5^) one 
has (neglecting again V{if)) 



Ssi 



diifdjtp 



XiXj 1' dtp 
du 



With h. 



finds 



p, - -T^P^ 

and the gravitational energy density Pg reads 



(62) 



(63) 



^''For a single cosmon lump the sign Al/ may be positive or 
negative, depending on the sign of dtp/du. 



^^Alternatively, one may translate the metric 
monic coordinates. 



into har- 



9 



> 



(64) 



This is to be compared with the energy density in the 
cosmon field 



Pc 



1 



(65) 



The average in eqs. (|6J) and (|65|) is given as an integra- 
tion over the volume of the lump 



< T{u) 3 



duu^T{u)/(p\j - ul 



(66) 



where Uh corresponds to the radius of the "bulk" of 
the galaxy and must be larger than the critical value 
Uc = Rn/e for the central singularity. (We recall here 
that the averaging needs to be done with respect to the 
background metric g^i/ such that the volume is just the 
cartesian volume in the coordinates x. It does not in- 
volve the "microscopic volume" which would have an ad- 
ditional factor ^ = ^1/2(73/2.) 

We note that C{u) becomes smaller than one inside 
the halo such that pg is indeed negative and Pg positive. 
(The numerical prefactors will be altered if eqs. (|2l| ) 
and ( pO| ) are used instead of eqs. ( p3| ) and (26).) This 
demonstrates how the cancellation between positive pc 
and negative pg becomes visible already in the linear ap- 
proximation. With i? w 1 and 



dip 
du 



4M2 



In [u/uc) 



(67) 



we find that the integrands relevant for pg and pc^ re- 
spectively, can be characterized by 



dip 



4M2 J 



u In (u/itc) 



dip 
du 



2, ,3 



Rlln^u/u,) 



(68) 



They are both dominated by the region u Uf,. We 
conclude that the energy density and pressure are actu- 
ally dominated by the interior of the halo and/or by the 
bulk. An assumption about a halo domination is actu- 
ally not justified, nor is the linear approximation for the 
computation of pg and pg . 

Nevertheless, the need of a large cancellation between 
a positive cosmon energy denstiy and a negative gravi- 
tational energy density remains true for large I7I, irre- 
spective of the shortcomings of the above calculation. 
Already the integration of the cosmon energy density 
over the halo exceeds the total mass by a large factor 
> Rh /Rs ~ I7I ■ The total sum ( ^ ) can only be balanced 
by a negative gravitational energy density of almost equal 
(averaged) size! We may summarize our discussion be ex- 
tracting the following general features for large cosmon 



fluctuations: The cosmon energy density pc is positive 
and the pressure Pc negative, typically with pc ~ —pc/3. 
This is accompanied by a negative gravitational energy 
density pg and positive gravitational pressure Pg. For 
static isotropic configurations large cancellations occur 
both for Pc + Pg and pc+Pg, implying w -Pg/S. (This 
differs from eq. (^) which involves unjustified approx- 
imations.) For more general, in particular non-static, 
large cosmon fluctuations the detailed balance between 
gravitational and cosmon contributions may not occur 
anymore. It is plausible, however, that the above flnd- 
ings about the sign of the various contributions remains 
valid. 

A very simple, but perhaps important, observation 
states that the cosmon pressure pc is likely to be neg- 
ative. Indeed, for large fluctuations we may neglect the 
subtraction of the potential and kinetic energy of the 
background field ip. The cosmon pressure is then given 

by 

p, ^V{ip) - ^9''d,^d,ip + > (69) 

We observe a negative contribution from the gradient 
term reflecting spatial inhomogeneities of ip. Also the 
contribution of the potential is negative and only a fast 
time variation could cancel these two negative contribu- 
tions. 

Imagine now a period in the cosmological evolution 
where the cosmon fluctuations become substantial and 
their negative pressure is not (or only partially) cancelled 
by the pressure of metric fluctuations. The cosmologi- 
cal evolution would then be substantially affected by the 
negative pressure of cosmon dark matter. Furthermore, 
the cancellation between cosmon and gravitational en- 
ergy density could be more effective than for the pressure. 
This could lead to a situation where the total energy 
momentum tensor is dominated by cosmon dark matter 
and quintessence with a substantially negative equation 
of state w 



= P/P 



Pc 



■Pg+Ph 



Pc + Pg + Ph 



(70) 



In fact, the pressure of dark energy ph could also turn 
negative if the potential dominates over the kinetic en- 
ergy during such an epoch. If w becomes smaller than 
-1/3, the expansion of the universe accelerates 



a 



P 



12 AP 



(1 + 3w) 



(71) 



It is tempting to speculate that such a situation may oc- 
cur towards the end of structure formation. The dom- 
inant contribution to Pc + Pg may arise from cosmon 
inhomogeneities on the scales of clusters or larger. It 
is even conceivable that the present acceleration occurs 
only effectively for the metric averaged over a volume 
corresponding to a redshift z of the order one. At earlier 
terms it may have been "visible" in the averaged metric 
relating to a smaller effective volume. 
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VII. CONCLUSION 



As a conclusion, let us turn back to the question asked 
in the title: can structure formation influence the cos- 
mological evolution? We have presented in this paper a 
few estimates and simple model calculations within a for- 
malism which describes the backreaction of fluctuations. 
We find it unlikely that standard cold dark matter fluc- 
tuations lead to a substantial effect, even though these 
fluctuations are today large and strongly nonlinear. The 
basic reason is that the averaged Einstein equations are 
linear in the energy momentum tensor. The direct contri- 
bution of fluctuations in the energy density and pressure 
therefore cancels by virtue of the averaging. An indirect 
effect of these fluctuations shows up in the form of in- 
duced metric fluctuations. This effect is related to the 
gravitational energy density and pressure. We have seen, 
however, that the size of this induced fluctuation effect 
is small unless a substantial part of the matter is in re- 
gions with strong and time-varying gravitational fields. 
Furthermore, we have seen that by a "cosmic virial the- 
orem" the gravitational pressure cancels the effect of the 
pressure of cold dark matter. 

The situation can change drastically in presence of a 
scalar "cosmon" field mediating quintessence. If the cos- 
mon fluctuations grow large, their contribution to the 
backreaction becomes typically quite sizeable. The av- 
eraged Einstein equation, as well as the averaged scalar 
evolution equation, are not linear in the cosmon fluctu- 
ations. In contrast to standard cold dark matter large 
fluctuations make therefore directly a large contribution 
to the averaged equations. Our computation for cosmon 
lumps has revealed that typically the induced gravita- 
tional energy density and pressure are also large. For 
a collection of static and isotropic cosmon lumps this 
"gravitational backreaction" cancels the "cosmon back- 
reaction" to a high degree. For the pressure one observes 
a matching of a negative cosmon and a positive grav- 
itational contribution. For more general large cosmon 
fluctuations, in particular if they are not static, this can- 
cellation may not be perfect. A large backreaction effect 
would then be expected for large cosmon fluctuations. 
We conclude that the backreaction could substantially 
influence the cosmological evolution after the time when 
large cosmon fluctuations have developed. 

We have also argued that the equation of state of the 
combined cosmon and gravitational fluctuations may be 
substantially negative. In this event a growth of fluctu- 
ations in the cosmon fleld towards the end of structure 
formation could trigger an acceleration of the expansion 
of the universe and provide an answer to the question 
why such an acceleration happens "just now" . Many 
pieces of the scenario outlined here are, however, fairly 
speculative. In particular, it remains to be seen if a re- 
alistic effective action for the cosmon field can be found 
such that the cosmon fluctuations indeed grow large in 
consistency with present observational information. 
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